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Abstract. Based on extended multicolour CCD photometry of the triple-mode radial pulsator V823 Cas we studied 
the properties of the coupling frequencies invoked by nonlinear processes. Our results support that a resonance 
connection as suggested by Antonello & Aikawa (1998) affects the mode coupling behaviour. The Pl/PO period 
ratio of V823 Cas has an "out of range" value if compared with the period ratios of the known double mode 
pulsators, while the P2/P1 period ratio is normal. The periods and period ratios cannot be consistently interpret 
without conflict with pulsation and/or evolution models. We attempt to interpret this failure by the suggestion 
that at present, the periods of V823 Cas are in a transient, resonance affected state, thus do not reflect the true 
parameters of the object. The anomalous period change behaviour of the fundamental and second overtone modes 
supports this idea. We have also raised the possibility that a /o + /2 = 2/i resonance may act in triple mode 
pulsators. 

Key words. Stars: individual: V823 Cas - Stars: variables: Cepheids - Stars: oscillations - Stars: evolution - 
Techniques: photometric 



1. Introduction 

Stellar pulsation often manifests multiperiodic behaviour. 
The simultaneous excitation of different modes is a com- 
mon property of e.g. S Set variables. In these stars both ra- 
dial and nonradial pulsation modes can be excited, and es- 
pecially in case of a dense frequency spectrum with many 
nonradial modes, the identification of the modes still has 
not been settled. 

Multiperiodicity of radial mode pulsators is also ob- 
served in different types of variables. Radial mode pulsat- 
ing variables are important targets of astero-seizmology 
as mode identification is unambiguous in many cases. 
Consequently, the pulsation properties of these stars re- 
veal the physical parameters quite precisely by using pul- 
sation models and/or empirical relations developed for the 
different types of variables along the instability strip. Two 

* Tables 1 — 5 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
or via_http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/_ 



radial modes can be excited simultaneously in Cepheids, 
HADS (high amplitude 6 Set stars) SX Phe and RR Lyrae 
variables as well. These, so called double mode pulsators, 
oscillate either in the fundamental and first overtone 
modes or in the first and second overtones. 

In three stars; AC And, V823 Cas and V829 Aql, 
the three lowest radial modes have been detected to 
be excited simultaneously. Very recently two triple- 
mode shor t period Cepheids were also discovered in 
the LMC l)Moskalik fc Dziembowskil l2005(l . these stars 
oscillate, however, in the first three overtone modes. 
As a result of the mass photometry of t he LMC and 
SMC ( OGLE, MACHO EROS p roiects. lUdalski et alJ 
( I999I): ISoszvnski et all l|200(]() : lAlcock et alJ l|l999|) : 
iBeaulieu et al.l l)l997|) ') dozens of short period (P < 1 d) 
double mode Cepheids were discovered, both fundamental 
and first overtone {FU / FO) and first and second overtone 
(FO/SO) pulsators. These results have also shown that 
in the LMC and SMC most of the P < 1 d Cepheids have 
double mode properties, indicating that in their physi- 
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cal parameter regimes multi-mode pulsation is favoured. 
Because of the large distance and the consequent faint- 
ness of the LMC/SMC multi-mode pulsators, their inves- 
tigations are strongly restricted. The three galactic triple- 
mode radial pulsators are also sparsely examined, in spite 
of that our knowledge on stellar pulsation would highly 
benefit from their thorough observational and/or theoret- 
ical studies. 

According to their fundamental periods AC And and 
V823 Cas {Pq ~ 0.7 d) resemble to metal poor RR Lyrae 
stars, while V829 Aql with Pq = .29 d is classif i ed as a 
post main sequence S Set star bv iHandler et alJ l)l998() . 
The metal content of AC A nd is, however , extremely high 
for a t ypical RR Lyrae star llPrestonll 9*591) . iFitch &: SzeidH 
l)l97 6'l and 'Kov acs fc BuchleJ l)l994j) based on pulsation 
and evolution models have concluded that it is most prob- 
ably a higher mass evolved star similar to S Set variables. 
iFernid l)l994|) suggested that AC And may be the missing 
link between S Set stars and classical Cepheids. 

V823 Cas was disco vered to be a. triple mode radial 
pulsator by S. Antipin (l Antipin|ll997f) using the Moscow 
photographic plate collection taken with the 40cm astro- 
graph in Crimea between 1948 and 1995. Data showed 
without doubt light variation with three independent pe- 
riodicities matching most probably to the three lowest ra- 
dial mode oscillations. Until today, no further information 
on this important object has emerged, therefore, we de- 
cided to carry out an extensive multicolour photometric 
investigation of V823 Cas. 

2. Observations, data reduction 

The observations were obtained with the automated 
60cm telescope of the Konkoly Observatory (Svabhegy, 
Budapest) equipped with a Wright 750x1100 CC D (pa- 
rameters and calibration are given in iBakosI ljl999l) ) using 
BV{RI)c filters. The field of view was 17'x24'. Reduction 
processes were performed using standard IRAF^ pack- 
ages. Data were corrected for atmospheric extinction. 

About 3 600 frames were obtained in each passband 
on 38 nights between 25 September and 14 December in 
2003 (JD 2 452 908 - 2 452 988). Transformation to the 
standard system was done using the B,V, Rc, Ic magni- 
tudes of surrounding stars observed by A. Henden with the 
USNO Flagstaff Station 1.0 m telescope equipped with a 
SITe/Tektronix 1024x1024 CCD. A complete table of the 
positions and BV{RI)c magnitudes of stars in the field 
centred on V823 Cas, is given in Tabled (available elec- 
tronically from the CDS) . The following formulae were de- 
rived to transform the instrumental data to the standard 
system: 

V = 0.989W -I- 0.095(6 -w) 

B = 0.987&- 0.032(6 -w) 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



Table 1. electronic table: BVRI magnitudes of calibrating 
stars in the field of V823 Cas 

Table 2. electronic table: V823 Cas B data 

Table 3. electronic table: V823 Cas V, B~V, V~Rc, V-h 
data 

Table 4. electronic table: V823 Cas R data 



Rc = 0.986r + 0.112(w-r) 
Ic = 0.989i - 0.012(w - i) 

Aperture photometry was applied to measure the rel- 
ative magnitudes of V823 Cas to CSC 04018-01777 {V = 
ir."235, S-F = l?n73, V-Rc = 0^617, V-I^ = 1^202). 
HD 134, CSC 04018-01891 and CSC 04018-01661 were 
used as check stars, neither of them showed any signifi- 
cant variability relative to GSC 04018-01777. 

Photometric data are available elec- 
t ronically at the CD S. In Table El - El 
( http: / / cdsweb.u-strasbg.fr /cgi-bin/qcat? J / A-|- A/ ) 
Column 1 lists the HJD of the observations, and 
Column 2 gives the differential magnitudes of V823 Cas 
with respect to GSC 04018-01777 for the B,V,RcJc 
colours. B — V, V — Rc, and V — Ic colour curves were 
also derived utilizing the V observations and fitted values 
of the B, Rc, Ic curves according the Fourier solutions 
given in Table jSl and Table [7\ for the moments of the V 
measurements. The colour indices are given in the 3. ,4., 
and 5. Columns of Tabled 

3. Photometric results 

3.1. Frequency components, amplitudes, and phases 

The light curve solution given in Table El and Table El 
was determined from the Fourier analysis of the photo- 
metric data using the utilities of the program package 
MUFRAN (KoUath 1990). Nonlinear regression facilities 
of Mathematica (Wolfram Research, Inc) were also ap- 
plied. The listed frequencies fit the B,V, Rc, Ic data with 
0.0126, 0.0079, 0.0081, 0.0079 mag r.m.s. scatter, respec- 
tively, which is about the observational accuracy limit. 
The V light curve and the fit are shown in Fig. ^ 

Besides the frequencies of the three dominant radial 
modes and their harmonics, 26 of the possible linear com- 
bination coupling frequencies could be identified in the 
successively prewhitened spectra of the B data. 

Frequency component identification was performed on 
the B dataset as the amplitudes of the signals are the 
largest in B. In order to refine the /o,/i,/2 frequency 
values least squares solutions were calculated within the 

Table 5. electronic table: V823 Cas / data 
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Identification 




frequency 


Ab 


Av 


Ar^ 


Ai^ 


Ab-v 


Av-Ra 


Av-i, 








cycle/day 








mmag 










/i 




1.950515 


176.45 


130.23 


102.37 


79.59 


46.42 


28.01 


51.22 




2/i 




3.901030 


37.15 


28.19 


21.97 


17.17 


8.77 


6.31 


11.09 




3/i 




5.851545 


13.13 


10.25 


7.67 


6.96 


3.00 


2.54 


3.25 




4/i 




7.802060 


3.62 


3.08 


2.52 


2.27 


0.50 


0.60 


1.11 


fo 






1.494765 


115.95 


86.00 


67.10 


52.28 


30.17 


19.26 


34.60 


2/(1 






2.989530 


16.64 


12.28 


9.78 


8.09 


4.41 


2.49 


4.20 


3/o 






4.484295 


2.14 


1.30 


0.88 


0.99 


1.31 


1.29 


1.31 






f2 


2.433590 


29.16 


22.04 


17.67 


13.57 


7.21 


4.37 


8.47 


fo 


/i 




3.445280 


44.77 


34.02 


26.97 


20.54 


11.23 


6.87 


13.22 


-fo 


/i 




0.455750 


25.10 


19.45 


15.45 


12.15 


5.67 


4.02 


7.28 


2/o 


/i 




4.940045 


11.76 


8.79 


7.30 


5.99 


3.17 


1.47 


2.74 


2/o 


-/i 




1.039015 


5.26 


3.72 


3.74 


2.68 


1.70 


0.44 


1.01 


3/o 


/i 




6.434810 


3.44 


2.12 


1.56 


1.56 


1.28 


0.62 


0.77 




/i 


f2 


4.384105 


10.11 


8.07 


6.31 


4.43 


2.29 


1.79 


3.73 




-/i 




0.483075 


6.94 


5.43 


3.88 


3.35 


1.43 


1.64 


2.13 


fo 


/i 


f2 


5.878870 


4.87 


3.18 


1.93 


2.39 


1.69 


1.40 


0.96 


-fo 


/i 


f2 


2.889340 


4.16 


3.21 


2.73 


2.06 


0.93 


0.63 


1.21 


fo 


-/i 


f2 


1.977840 


4.16 


3.24 


2.52 


2.15 


1.19 


0.94 


1.55 


2/o 


/i 


f2 


7.373635 


1.56 


1.52 


1.33 


0.76 


0.25 


0.27 


0.80 


fo 


2/i 




5.395795 


10.70 


7.91 


5.97 


4.51 


2.78 


2.10 


3.41 


-fo 


2/i 




2.406265 


13.50 


10.48 


8.67 


6.56 


2.97 


1.88 


4.00 


2/o 


2/i 




6.890560 


3.82 


3.05 


2.27 


1.89 


0.68 


1.03 


1.31 


-2/o 


2/1 




0.911500 


2.07 


1.92 


1.25 


0.89 


0.17 


0.70 


1.24 




2/i 


/2 


6.334620 


5.89 


4.18 


3.22 


2.61 


1.71 


0.97 


1.63 


fo 


2/i 


/2 


7.829385 


2.52 


2.09 


1.56 


1.46 


0.42 


0.55 


0.95 


-fo 


2/i 


/2 


4.839855 


2.11 


1.14 


0.73 


0.22 


1.31 


0.33 


0.76 


fo 


3/i 




7.346310 


6.33 


4.80 


3.56 


3.20 


1.65 


1.19 


1.64 


-fo 


3/i 




4.356780 


3.50 


3.15 


2.63 


1.55 


0.39 


0.65 


1.76 


2/o 


3/i 




8.841075 


1.72 


1.58 


0.63 


0.53 


0.52 


0.95 


1.15 




3/i 


/2 


8.285135 


3.05 


2.19 


1.27 


1.62 


1.00 


0.91 


0.59 




3/i - 


■/2 


3.417955 


2.16 


1.57 


1.56 


1.29 


0.50 


0.19 


0.31 


fo 


4/i 




9.296825 


2.89 


1.82 


1.65 


1.25 


1.19 


0.29 


0.58 


fo 




/2 


3.928355 


6.01 


4.22 


3.17 


3.36 


1.88 


1.11 


0.80 


-fo 




/2 


0.938825 


2.45 


2.22 


1.12 


1.16 


0.71 


1.12 


1.47 


/3 






0.817000 


4.02 


2.27 


1.87 


1.43 


1.73 


0.49 


0.89 



vicinity of the suspected values of /o, /i, /2, their harmon- 
ics and the detected hnear combination term frequencies. 
Linear combination frequencies were always set according 
to the actual values of the main frequency components. 
Solution with the smallest residual scatter was accepted 
for the fa, fi, /2 frequency values. The 2a asymptotic er- 
ror estimates of the corresponding periods are 0.000045, 
0.000013 and 0.000060 days, respectively. 

The formal errors (Icr) of the amplitudes are about 
0.2 — 0.4 mmag, the errors of the phases of signals with 
amplitude larger than 1 mmag are reliable within 0.1 — 10°. 

All but 3/i — 2/o 2-term coupling frequencies involving 
/i,2/i,3/i, and /o, 2/o appear in the spectra, and six 3- 
term linear combination frequencies are also detected. The 
amplitudes of the coupling components are always smaller 
than that of its compounding frequencies with one excep- 
tion. The amplitude of the /i -I- 3/o component is higher 
than the amplitude of the 3/o in each colour. Actually 



the 3/o harmonic component is only marginally present 
in the data, we involved this component to the solution 
only because a real frequency component was found at 
6.43481 c/d, which we could interpret only as /i -I- 3/o. 
However, considering the marginal amplitude of the 3/o 
component, this interpretation might not be correct. 

After the removal of the radial modes, their harmonics 
and the coupling terms from the data, the residual spectra 
in each colour show evidence of a component at 0.817 c/d 
as shown in Fig. El This frequency, denoted as /a in TableEl 
and Tabled cannot be resolved as the linear combination 
of the radial modes. As no indication of this frequency was 
found in the comparison - check stars data, we suppose it 
as a real frequency component of V823 Cas. Taking into 
account the small amplitude and the relatively long, 1.22 d 
periodicity of this variation one may speculate about its 
origin as a high-order gravity mode, similarly as in the 
7 Dor variables. In order to check the reality of this ex- 
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Table 7. Phases of the detected frequencies in B ,V, Rc, Ic, B — V,V — Rc,V — Ic using cosine term Fourier sums and 
To = 2 452 908 initial epoch value. 



Identification 




frequency 






















cycle/day 








deg 








— 


fi 






—85.24 


—86.21 


—87.87 


—90.72 


—82.60 


—80.09 


—79.14 




2/i 






17.17 


16.42 


15.56 


14.11 


20.43 


18.79 


19.82 


— 


3/i 






167.96 


165.83 


165.24 


167.21 


179.39 


166.14 


161.51 


— 


4/i 




1 . U w 


—40.83 


—40.81 


—36.70 


— 24.62 


—46.76 


—57.82 


— 74.12 


/o 


— 






—9.53 


— 11.77 


— 14.50 


— 18.37 


—3.11 


—2.26 


— 1.81 


2/o 


— 






220.48 


218.75 


217.74 


215.09 


225 52 


222.46 


225.57 


3/o 


— 






69.47 


66.47 


57.45 


58.57 


247.29 


66.61 


66.60 


— 


— 


.12 


ij . t: tj ij i,/ 1^ 


—32.64 


—33.35 


—34.37 


—34.62 


—30.66 


—29.49 


—31.23 


fo 


/i 




44^1280 


156.53 


156.86 


157.49 


157.50 


156.46 


154.16 


155.79 


-/o 


/i 




L/ .^iijij ( iJ\J 


124.69 


124.51 


123.20 


125.22 


124.11 


130.23 


123.71 


2/o 


/i 






16.59 


18.51 


20.08 


21.15 


7.98 


14.15 


14.22 


2/o 


-/i 






124.89 


131.37 


138.78 


130.44 


108.37 


43.21 


135.37 


3/o 


/i 






259.64 


248.62 


244.81 


234.59 


275.08 


259 32 


277.82 


— 


/i 


fo 




78.22 


81.44 


82.74 


80.75 


55.23 


82.75 


85.81 




-fi 


./2 




234.69 


235.40 


230 99 


238 92 


229.87 


246.51 


230 05 


fo 


fi 


./2 


5 87SS7n 


—78.41 


—78.36 


—62.65 


—65.54 


—79.30 


—98.87 


— 108.18 


-fo 


fi 


^2 


2 88Q'^4n 


255.21 


262 08 


253 65 


252 93 


235 83 


295 82 


273.64 


fo 


-fi 


fo 


1 q77S4n 


54.17 


67.07 


57.34 


46.21 


10.08 


96.75 


98.38 


2/o 


fi 


fo 
j2 


7 S7'Hfi'H.'i 


119.60 


131.29 


137.83 


145.01 


44.71 


91.55 


116.02 


fo 


2/i 






—80.14 


—77.50 


—70.74 


— 72.81 


—88.61 


—96.49 


—83.44 


-fo 


2/i 




2 4nfi2fifi 


242.97 


240.53 


238 93 


237.65 


249.56 


249.04 


245.89 


2/o 


2/i 




fi 8Qn'ififl 


158.83 


158.21 


171.31 


170.19 


169.90 


126.06 


138.32 


-2/o 


2/i 




0.911500 


-3.69 


2.88 


3.42 


-27.58 


-77.19 


0.01 


22.05 




2/i 


/2 


6.334620 


224.63 


224.04 


228.02 


217.36 


225.35 


209.85 


234.81 


fo 


2/i 


/2 


7.829385 


95.76 


91.42 


94.31 


114.91 


112.82 


84.46 


55.10 


-fo 


2/i 


/2 


4.839855 


72.13 


78.17 


67.14 


108.60 


71.59 


104.80 


66.94 


fo 


3/i 




7.346310 


43.06 


41.77 


41.87 


49.15 


43.61 


42.78 


28.25 


-fo 


3/i 




4.356780 


30.46 


39.13 


33.64 


38.47 


-44.35 


54.05 


35.88 


2/o 


3/i 




8.841075 


-81.99 


-67.51 


-65.18 


-37.29 


-147.19 


-68.90 


-79.12 




3/i 


/2 


8.285135 


29.87 


20.70 


18.68 


21.75 


54.73 


21.53 


13.96 




3/i - 


-/2 


3.417955 


34.45 


48.91 


52.36 


43.12 


-1.11 


-12.86 


58.26 


fo 


4/i 




9.296825 


199.65 


185.55 


178.67 


191.70 


218.01 


236.97 


175.97 


fo 




/2 


3.928355 


134.81 


135.73 


141.33 


137.49 


140.06 


114.45 


123.12 


-fo 




/2 


0.938825 


171.24 


155.66 


169.61 


197.53 


222.10 


141.50 


124.41 


/3 






0.817000 


76.96 


82.20 


75.96 


76.14 


68.92 


108.98 


92.31 



planation, however, further observational and theoretical 
confirmations are needed. 

No other frequency in the residual spectrum was de- 
tected, consequently the solution given in Table El and 
Table [T] is a full description of the data within the lim- 
its of the accuracy of the observations. 

In most cases the linear combination components in- 
volving negative term have smaller amplitudes than the 
combination component of the same but positive terms, 
in agreement with model results, that nonlinear pulsa- 
tion favours combi nation frequency terrns wh ere both i 
and j are positive ijAntonello fc AikawalllQQ^ . There is 
one significant counter-example, 2/i — /g has larger am- 
plitude than 2fi + fo in each colour. The amplitudes of the 
3-term combination frequencies of the same components 
have similar amplitudes independently of the signs. 



The amplitude ratios of the combination frequencies 
defined as 

R^Jk = + \j\ + \k\)A,jk/{\i\Awo + \Mow + \k\Aoai) 
decrease exponentially with increasing order (|i| + |i| + |fc|) 
as shown in FigO Aijk denotes the amplitude of the ifo + 
jfi + kf2 frequency component, i.e., Aiqq = A{fo), Aqiq = 
and Aooi - A(/2). 

The amplitudes of all the detected frequencies decrease 
towards longer wavelengths. Most of the frequencies of the 
light curve solution appear in the colour indices as well, 
typically with the largest amplitude in the V — Ic colour. 

The phases of the three pulsation modes tend to be 
smaller towards the longer wavelengths. The phase differ- 
ences between the B and Ic data for fo, fi and /2 are 
9°, 5°, and 2°, respectively, while the phases of the colour 
indices are typically larger than the phases of the light 
curves, e.g. the phases of /o, /i and /2 of the B — V curve 
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Fig. 1. V light curve of V823 Cas fitted with the Fourier 
solution given in Table El and Tabled 



are larger by about 8.7°, 3.6°, and 2.7°, than the phases of 
the V light curve solution. In contrast with these system- 
atic trends of the phases of the three radial modes, the 
phases of the two largest amplitude linear combination 
terms, /i -|- /o and fi — fa, are identical in each wave- 
length bands. The phases of the smaller amplitude linear 
combination terms are not accurate enough to detect any 
re al wavelenght dep e ndent trend in the phase behaviour. 

IPoretti fc Pardol (^23) have shown that in galac- 
tic double-mode Cepheids the different order generalized 
phase differences (G^fe, for a definition see .Antonellg 
l|1994() ) of the linear combination term frequencies tend 
to have similar values. In Table |S1 the different order gen- 
eralized phase differences of our B light curve solution and 
also the Rijk amplitude ratios are listed, in the order of 
decreasing amplitude. The la errors of the phases are also 
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Fig. 2. Residual spectra of the B,V,Rc,Ic data after 
prewhitening with the frequencies of the three radial 
modes, their harmonics and the linear combinations of 
these frequencies appearing in the Fourier spectra. In each 
colour the highest signal appears at 0.817 c/d indicating 
the detection of a real frequency component. 



Table 8. Generalized phase differences and amplitude ra- 
tios of the frequencies of the B light curve solution: 
Gijk = 4>tjk - {i4'im + j0oio + ^0001 ), 
R^Jk = + lil + \k\)A,jk/{V\Awa + IjI^oio + |fc|^ooi). 
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given, in order to easily estimate the reliability of the Gijk 
phase difference values. For compatibility with other sim- 
ilar data, the phase differences are given in radians. The 
averages of the Gijk phase differences are 3.55, 1.07, 5.01 
(5.26 omitting the most deviant 2.577 value), and 3.34 for 
the 2., 3., 4. and 5. order terms, respectively. These values 
are systematically smaller by about 0.5 — 1.0 rad as the 
corresponding 4.31, 2.17, 6.24, 3.85 averages of the Gij 
values of galactic double- mode Cepheids calculated from 
the data listed bv lPoretti fc Pardol l)l997l) . Taking into ac- 
count the decreasing trend of the phase differences with 
decrea sing periods as shown in Fig. 2 in lPoretti fc Pardol 
l|l997j) the short periods of V823 Cas may account for its 
smaller Gijk values. 

As indicated by lAntonello fc Aikawal l)l998|) the only 
plausible explanation that can change the phases and can 
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Fig. 3. Amplitude ratios of the different order combina- 
tion frequency terms. Filled and open circles denote fre- 
quencies with positive i, j, k values and involving negative 
terms as well, respectively. The decrease in the amplitude 
ratios is exponential with increasing order. 

lead to coherent values of the different order phase dif- 
ferences is a resonance between the pulsation modes. In 
Fig 0] the B light curve is folded according to the periods 
of the three radial modes (top panels), after prewhiten- 
ing with the frequencies and harmonics of the other two 
modes (middle panels), and the residual light curve af- 
ter the removal of all the three radial modes' light curves 
(bottom panels). These plots still show some kind of reg- 
ular behaviour according to the main periodicities. It is 
also evident that extreme large amplitudes occur when 
each of the modes are simultaneously around maximum 
phase indicating also that some kind of resonance interac- 
tion may be responsible for the occurrence of the coupling 
frequencies. 

4. Comparison with model results 

Period ratios of different type radial multi-mode pulsators 
has been collected and the result is summarized in FigO 
The ~ 0.01 separations of the Pi/Po sequences of galactic, 
LMC and SMC Cepheids reflect the different metal con- 
tent of the stars {Z = 0.02, 0.008 and 0.004, respectively). 
The Pi /Pq ratio of V823 Cas is evidently discrepant from 
the global log Pi ~ Pi/ Pq relation defined by the different 
types of variables. This is most probable because of an 
anomalous value of Pq, as the P2/P1 period ratio fits well 
the observed values of other multi-mode pulsators. 

In order to find an answer for the anomalous 
value of the fundamental mode period first we have 
checked whether the observed periods and period ra- 
tios can match stellar pulsation and evolution model 
results. The difference between pulsation periods de- 
rived from linear models and the real nonlinear periods 
are supposed to be negligible. Thus linear periods 
are widely used to determin e physical parameters of 
multi- mod e pulsators fe.g. iMoskalik k, Dziembowskil 

i^ili IPetersen Christensen-DalsgaaTd! ill flflfJl : 

lAntonello fc Aikawal |l£2a))- iKovacs fc Biichleil l|l994h 
(hereafter KB94) calculated an extended set of linear 
radial models to explain the periods of AC And, the 
other triple-mode variable with very similar periods. We 
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Fig. 4. B light curve of V823 Cas folded with the peri- 
ods of the three radial modes (top panels). The synthetic 
light curves corresponding to the three periods are also 
shown calculated from the amplitudes and phases of the 
three radial modes and their harmonics according to the 
data given in Table El and Table |3 Middle panels show 
the folded light curves after subtracting the other two ra- 
dial modes and their harmonics. The shape of the light 
curve of a given mode is strongly variable depending on 
the relative phases of the modes. In the bottom panels 
residual data after removing the light curves of the three 
radial modes are shown. These panels clearly indicate that 
the amplitude of the light variation increases significantly 
when each of the modes are simultaneously around max- 
imum brightness. This enhancement may account for the 
numerous coupling frequencies appearing in the Fourier 
spectru m of the data in accordaii ce with the model re- 
sults of lAntonello fc Aikawal l(l998^ . 



have compared the pulsation periods calculated from 
the KB94 formulae with those given in the evolutionary 
and pulsation models of ijAlibert et alJll999|) . The period 
of a given model agrees within 1% accuracy with the 
period calculated from KB94 models (within the common 
parameter regimes) indicating consistency of the two 
model groups. 

The observed periods, and period ratios Po — 
0.669001 d, Pi/Po = 0.7663, P2/P1 = 0.8015, can be fit- 
ted with period ratio distance A < 0.003 (for definition see 
KB94) accuracy with certain high-mass-sequence 2 mod- 
els of KB94. Solutions can be found for the Z=0.004 and 
Z=0.01 models with 3.7 - 4.0 Mq and 3.85 - 4.25 Mq 
mass, respectively. The metal-poor solutions [Z = 0.004) 
are at too high, T > 7 000 K temperature out of the in- 
stability strip, with negative growth rates of the detected 
modes, therefore we exclude these among the real possi- 
bilities. There is also indication that solutions might also 
exist for solar composition, Z=0.02 models, but at 'out 
of range' M ~ 4.5 mass value, which was not covered 
with the KB94 models, thus these pulsation solutions are 
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Fig. 5. Period ratios as a function of log Pi of different 
types of multi-mode radial pulsators. Though Pi is ac- 
tually the same for AC And and V823 Cas, their Pi/Pq 
period ratios differ about 0.03 which cannot be explained 
with any major difference in their physical parameters 
consistently with pulsation and evolutionary models. 

somewhat uncertain. Moreover, considering the Period - 
T emperature connec t ions v alid for Ceph eids, e.g. , Fig. 1 
in lBaraffe fc AlibertI l|200l[l and Fig. 3 in lSandagd l|2004|) . 
the temperature of Cepheids with periods shorter than 
1 day must be higher than 6000 K. All the Z=0.02 model 
period solutions fall below 6000 K, thus these are not valid 
solutions as well. 

The parameter regimes of the possible solutions are 
shown in Fig. The physical parameters of these models 
seems to be, however, in conflict with evolutionary results. 
At a given mass and chemical composition any canonical 
evolutionary model predicts the possible lowest luminosity 
within the instability strip during its 'first crossing', evolv- 
ing off the main sequence. However, each of the appropri- 
ate evolutionary rnodels llAlibert et al.lll99^ iBono et all 
l2000HGirardi et al.l2000l:lYi et alJ200^ cross the instabil- 
ity strip at significantly larger luminosities than the solu- 
tions of the periods indicate even during the first crossing. 
The minimal luminosity at a given metallicity and mass 
values according to evolutionary models are also indicated 
in Fig.El The 100 — 600 Lq discrepancy between luminosi- 
ties allowed by stellar evolution and luminosities derived 
from the pulsation solution of the periods indicates also 
that the periods of V823 Cas are anomalous. 

This problem arises from the unusually large value 
observed for the Pi/Pq ratio (0.766). At Pq w 0.7 d 
the normal value of Pi/Pq would be less than 0.75 
for any consistent model. The extended linear nonadi- 
abatic pulsati on models for galactic, LMC and SMC 
compositions l|Morgan fc WelchI Il997^ strengthen that 
Pi/Pq can never be larger than 0.75 in any real 
case. There are two ways to increase the period ra- 
tio at a given period, towards more metal poor mod- 
els, and/or with models of smaller L/M ratio the Pi/Pq 
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Fig. 6. Comparison of pulsation solutions for the observed 
periods and period ratios of V823 Cas with evolution- 
ary model results. Filled circles denote pulsation solu- 
tions according to the KB94 models for Z=:0.004 (light 
gray), Z=0.01 (black) and Z=0.02 (dark gray) composi- 
tions, respectively. Other symbols indicates the luminosi- 
ties of evolutionary tracks during first crossing of the same 
mass and metallicity models at the temperature of the pul- 
sation solutions. Second, third crossing luminosities are at 
even larger luminosities. The corresponding evolutionary 
models have the same mass within ±0.25 Mq and the 
same metallicity within Z — ±0.002 ranges as the pulsa- 
tion model solutions have. Each of the evolutionary models 
indicates larger luminosity than the pulsation solution. 

ratio is larger l| Pe1;;^r^e n fc Christensen-DalsgaardI Il99fil: 
iKovacs fc Buchleilll994|) . For V823 Cas. however, both of 
these possibilities can be excluded, as already mentioned 
the metal poor period solution fall out of the instability 
strip, and small L/M ratio is not allowed by evolutionary 
models. 

As a conclusion, we have found that the period ratios 
of V823 Cas cannot be explained consistently with stellar 
evolution and pulsation models. We suppose that this is 
not the consequence of any serious defect /inadequacy of 
the models but more probably arises from a temporal, 
transition behaviour of V823 Cas, which produces unusual 
period ratios (see the next section). 

5. Period changes 

The extended photographic data complemented with the 
recent CCD observations enable us to de rive relia.ble pe- 
riod changes of the three observed modes. lAntininl (^Q^) 
analysed two sets of the photographic data (between 
2 432 853-2 439 145 and 2 440 071-2 449 633, respectively) 
and detected period increase of the fundamental and first 
overtone modes. As photographic data are available al- 
ready from JD 2 429 076, we divided the photographic data 
into three subsets in order to get more details of the period 
change behaviour of V823 Cas. Table El lists the periods 



8 Jurcsik et al.: The triple-mode pulsating variable V823 Cas 



Table 9. Periods and period change rates of V823 Cas and AC And 



star JD-2 400 000 


^'o [d] 


2a{Po) Amp{Po) 


A [d] 


2a{Pi) Amp{Pi) 


P2 [d] 


2cr(P2) AmpiPi) 




V823 Cas 


52 908 - 52 988 


0.669001 


.000045 0.12 


0.512685 


.000013 0.18 


0.410916 


.000060 0.03 


CCD(B) (1) 


42 273 - 49 033 


0.668887 


.000005 0.13 


0.512660 


.000002 0.18 


0.411027 


.000005 0.05 


Pg (2) 


37 575 - 42 016 


0.668870 


.000009 0.12 


0.512646 


.000003 0.17 


0.411030 


.000007 0.05 


Pg (2) 


29 076 - 36 540 


0.668840 


.000009 0.13 


0.512628 


.000003 0.17 


0.411040 


.000007 0.07 


pg (2) 


AC And 


47 864 -49 039 


0.711253 


.000014 


0.525155 


.000009 


0.421088 


.000014 


Hip (3) 


36 459 - 37 949 


0.711232 


.000002 


0.525133 


.000001 


0.421072 


.000002 


pe(V) (4) 


30 258 - 35 009 


0.711227 


.000001 


0.525130 


.000001 


0.421067 


.000001 


pg (5) 


25 540 - 27 278 


0.711217 


.000011 


0.525108 


.000008 


0.421063 


.000016 


vis (6) 


Period change rates 


Po/Po 




Pi /Pi 




P2/P2 


[Myear^^] 





V823 Cas increasing 2.17 decreasing 

AC And 0.71 0.79 1.08 



References of the data: (1) this paper; (2) lAntipinI lll997f ): r3') lESAl ||1997D : r4') iFitch fc SzeidJ J1976D : f5 HGumanl lll982f) 76) 
FLausc (1932. 1933 193a) 



of the radial modes, their 2a error estimates and the pho- 
tographic or CCD(_B) amplitudes derived from nonlinear 
regression process of Mathematica. 

For comparison, the period change behaviour of 
AC And has been revised and listed in Table too. 
Periods for the different datasets of AC And has been 
determined similarly as described in Sect 3.1. 



I 2x10 'd 



2,430,000 



2,440,000 
JD 



2,450,000 



Fig. 7. Detected periods of V823 Cas at different epochs. 
Po - 0.66854 d (o), Fi - 0.51240 d (•), and P2 - 0.41088 d 
(A) are plotted for clarity. 



The period change rate of the dominant, first over- 
tone mode of V823 Cas and that of all the three modes 
of AC And well correspond to the theoretically expected 
evolutionary value during the first crossing of the insta- 



bility strip. The evolutionary period change rates during 
the second a nd third cro ssing are by two orders of magni- 
tude slower ijMoskalik fc Dzicmbowski 2005). The larger 
Pi/ Pi indicates somewhat larger mass of V823 Cas than 
AC And. 

In contrast, the periods of the fundamental mode 
and the second overtone of V823 Cas behave differently. 
According to the CCD observations the period change 
rates of these modes have been recently changed, Pq has 
increased, while P2 has decreased by about 0.0001 d dur- 
ing the last decades. The amplitude of the second overtone 
also shows significant changes, it has gradually decreased 
from 0.07 mag to 0.03 mag, in contrast with the ampli- 
tudes of the fundamental and first overtone modes, which 
differ only within the limits of the uncertainty for the dif- 
ferent sets of the observations. The reduction of the P2 
amplitude indicates that this mode may be just diminish- 
ing, and this is accomplished by the extreme period change 
behaviour. However, this effect cannot be solely account 
for the rapid period increase of Pq- As we have no reason 
to assume that canonical evolution and pulsation models 
which fail to interpret the period ratios of V823 Cas were 
in fault, we suggest that, instead, the fundamental mode 
period is somehow 'wrong'. 

An interesting possibility arises from the opposite sign 
period change of the second overtone mode. Extrapolating 
backward the linear period changes detected between 
JD 2 430 000-2 445 000, the period ratios should have been 
about Pi/Po = 0.76679 and P2/-P1 = 0.81089 about 3 300 
years ago. These period ratios fulfill the fo + h = 2/i res- 
onance criterion, as /o = 0.76679/i and = 1.23321/i, 
indicating that the three radial modes might have been in 
a resonance connection in a very recent past in astronom- 
ical time-scale. 

This resonance hypothesis would also help to under- 
stand the anomalous period ratios of V823 Cas. If we 
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assume that in the future the period change rate of the 
first overtone remains the same as detected for the last 50 
years, but the period increase of the fundamental keeps on 
to be larger than that as indicated by recent observations, 
i.e., 

A /-Pi - PolPo = C < 0, C - const. 
then, this would lead to change in the period ratio, as 

^^^^ - Pi/Po X (A/Pi - P^/Po) - C X Pi/Po. 
Consequently, with time elapsing, 

(Pi/Po)t = {Pi/Po)t=o X e^*. 
The exponential change in the period ratio results drastic 
changes during astronomically very short time base. The 
recent observations have shown that the period change 
rate of Pq has increased to 3 times that of the period 
change rate of Pi. Assuming the possible largest uncer- 
tainties, C is smaller than —2 Myear~^ at present, which, 
according to the above formula would lead to a decrease in 
the Pi/Po ratio from its unusual 0.766 to a normal 0.745 
value within 14 000 years. 

6. Absolute parameters derived from observed 
magnitudes, and P-L relations 

The magnitude and intensity mean V brightness and 
colours of V823 Cas using the standard magnitudes of the 
comparison star as given in Section 2 are the foUowings. 
Magnitude mean values: 

V = 11.175, B-V = 0.792, V-Rc = 0.473, V-Ic = 0.987. 
Intensity mean values: 

V = 11.169, B-V = 0.787, V~Rc = 0.471, V-Ic = 0.982. 

The interstel lar reddening in the d irection of V823 Cas 
according to the lSchlegel et al.l l(l998l) maps is quite large, 
E{B -V) = 0.887 mag. V823 Cas lies close to the galac- 
tic plane {b = 0.99°) to the direction of the Perseus Arm. 
The distance of the Perseus Arm is about 3 Kpc in the 
di rection of V82 3 Cas {I = 118°) as shown in Fig. 2 
in lOuillenl l(2002l) . Consequently, it can be assumed that 
within 3 Kpc interstellar reddening is much smaller than 
0.887 mag. Dereddened colours, using sta ndard redden- 
ing la w for BV{RI)c colours as given in ISchlegel et alJ 
1 1998^ fit the synth etic colours of atmosphere models 
I Castelli et all Il997j) according to the parameters hsted 
in Table Cni 

The observed periods of V823 Cas are too short for 
classical Cepheids and are too long for S Scuti models. 
This intermediate status makes comparison with model 
results uncertain, as in most of the cases conclusions can 
be drawn from extrapolation. Because of this deficiency 
and the uncertainties of data conversion between the ob- 
servable (magnitude, colour-indices) and the theoretical 
plane's (luminosity, temperature) quantities, the absolute 
parameters from the observed colours and periods cannot 
be determined accurately. 

A rough estimate can be, however, done by combining 
the possible solutions listed in Table [HH with predictions 
of empirical and theoretical Period - Luminosity relations. 
The 0.669 d fundamental mode, and 0.513 first overtone 



Table 10. Solutions for dereddening the observed 
coulours of V823 Cas. 



E{B-V) Z 


Teff 






Mv 








1.5 


2.0 


2.5 


3.0 [Kpc 


0.37 0.020 


6500 


-0.93 


-1.56 


-2.04 


-2.44 


0.19 0.010 


5750 


-0.35 


-0.98 


-1.46 


-1.86 


0.15 0.004 


5500 


-0.22 


-0.85 


-1.33 


-1.73 



periods correspond to My = —0.9 mag absolute visual 
magnitude during first crossing within ±0.3 mag range 
for Z = 0.02 — 0.004 matallicit y models according to the 
Per iod - Luminosity rela tions of iBaraffe fc AhbertI (|20^ 
and lAlibert et all ((l999l) . This value is in accordance with 
combined empiri cal P — L re lations of Cepheids and 5 
Scut i stars (e.g. fFerniel Il99d) Rec ent empirical studies 
e.g.. lUdalskil l)2000() : ISandagd ^^2004^ yield also consistent 
results within the given magnitude range. 

Comparing this value with the absolute magnitudes of 
different reddening and metallicity solutions at a distance 
of 1.5 — 3.0 Kpc as given in Tableliniwe conclude that if the 
Period - Luminosity relation of classical Cepheids is valid 
for V823 Cas, it is either a 1.5 Kpc distant, Tg// = 6500 K 
object with solar composition or a cooler, Z=0. 01-0. 004 
metallicity object at about 2 Kpc. As it was already men- 
tioned in Sect. 4, Cepheids with P < 1 day are hotter than 
6000 K, thus we conclude that the most probable solution 
for V823 Cas is E(B-V)=0.37, d=1.5 Kpc, T^fj = 6500 K, 
Mv = —0.9, Z = 0.02. Because of the uncertainties in- 
volved, this result has to be taken with caution. 

7. Summary 

The recent CCD observations of V823 Cas can be fitted 
with the periods of the 3 lowest radial modes, their har- 
monics and their linear combination terms within the ob- 
servational accuracy. No indication of nonradial p-modes 
has been found, but the probable existence of a small am- 
plitude signal with 1.22 d periodicity needs an explana- 
tion. 

We present, for the first time, extended multicolour 
observations of a multi-mode radial pulsator which is ac- 
curate enough to study the phase and amplitude relations 
of the radial modes and the different order coupling fre- 
quency components. These results give a deeper insight 
to the behaviour of the nonlinear processes taking place 
during stellar pulsation. 

The observed period ratios of V823 Cas cannot be con- 
sistently explained with canonical pulsation and evolution 
models. We interpret this anomaly as a consequence of a 
temporal behaviour of the periods which may originate 
from the very rapid evolutionary changes of the stellar 
parameters and/or a resonance connection between the 
radial modes' frequencies (/o + /2 = 2/i), which might 
have been fulfilled in the astronomically recent past. Both 
theoretical and observational investigations are needed in 
order to check the reality of this resonance during the 
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triple-mode phase of the pulsation, and to examine its ef- 
fect on the observed period values. 

If the actual values of the Pi and Pg period change 
rates remain unchanged, a normal Pi/Po period ratio is 
going to be reached in a short time, in about 10-20000 
years. This peculiar period change behaviour, together 
with the lapse of the P2 mode, as indicated by its ampli- 
tude decrease, may also hint that in its triple-mode state 
the periods of V823 Cas may have transient, resonance 
affected values, which do not reflect the real physical pa- 
rameters of the star. 

The unique period change behaviour of V823 Cas 
warns about the validity of drawing any conclusion about 
the evolutionary status and the direction of the evolution 
from the period change rates. A similar conclusion was 
reached from the unexpected period change behaviour of 
the two recently discovered over tone mode triple-mode ra- 
dial p ulsators in the LMC by iMoskalik fc Dziembowskil 
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